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It is well-known that the spacetime diagrams of some cel- 
lular automata have a fractal structure: for instance Pas- 
cal's triangle modulo 2 generates a Sierpinski triangle. Ex- 
plaining the fractal structure of the spacetime diagrams 
of cellular automata is a much explored topic, but virtu- 
ally all of the results revolve around a special class of au- 
tomata, whose typical features include irreversibility, an 
alphabet with a ring structure, a global evolution that is a 
ring homomorphism, and a property known as (weakly) p- 
Fermat. The class of automata that we study in this article 
has none of these properties. Their cell structure is weaker, 
as it does not come with a multiplication, and they are far 
from being p-Fermat, even weakly. However, they do pro- 
duce fractal spacetime diagrams, and we explain why and 
how. 



INTRODUCTION 

The fractal structure of cellular automata (CAs) has been a topic of in- 
terest for several decades. In many works on linear CAs, the authors 
present ways to calculate the fractal dimension or to predict the state 
of an arbitrary cell at an arbitrary time step, with much lower com- 
plexity than by running the CA step by step; However, their notions 
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of linearity are quite different. Often only CAs that use states in Z p are 
studied |*j Other approaches are more general, but still make certain 
assumptions on the time evolution or the underlying structure of the 
CA. In this work we try to loosen these restrictions as far as possible. 
We consider one-dimensional linear CAs whose alphabet is an abelian 
group. We show how they can be described by n x n matrices with 
polynomial entries and use this description to derive a recursion rela- 
tion for the iterations of the CA. This recursion relation enables us to 
formulate the evolution of the spacetime diagram as a matrix substi- 
tution system, which in turn gives us the means to calculate the fractal 
dimension of the spacetime diagram. 

Most of the methods we employ are commonly used in the study of 
CAs. To prove that the spacetime diagram converges we use elemen- 
tary graph theory and study the graph associated with the matrix sub- 
stitution system of the CA. The matrix substitution system is obtained 
by a recursion formula for the iterations of the CA and grouping of 
cells. 

Our interest in the fractal structure of CAs on abelian groups stems 
from our study of Clifford quantum cellular automata (CQCAs) ||T2| . 
We first noticed the self-similar structure while studying their long 
time behaviour [4, 3J. A CQCA maps Pauli matrices to tensor prod- 
ucts of Pauli matrices times a phase. If we neglect the phase, we can 
identify the Pauli matrices X, Y, Z with the elements of l\ via the 
mapping X i-S- (J), Y h-S> (J), Z h-» (°), 1 h-S> (°). Using this mapping 
we can simulate CQCAs with linear CAs on the alphabet ZS. The CAs 
corresponding to CQCAs have to fulfill rather strong conditions: they 
have to be reversible and preserve a symplectic form which encodes 
the commutation relations of the Pauli matrices |Tl2l . While our analy- 
sis is now much more general, our main example 0, whose spacetime 



diagram is shown in figure la is the classical counterpart of a CQCA. 

Our paper is organized as follows: in section [T] we give our defi- 
nition of a linear cellular automaton, introduce the formalism we will 
be working with, and state the main result: every linear cellular au- 
tomaton has a fractal structure. We will also introduce the example 
which will be the focal point of this article. In section|2j we will give an 
intuitive idea as to why the spacetime diagram of exhibits a fractal 



* We use the simple notation Z ( ; for the cyclic group of order d, instead of Z/ c?Z, as 
we are concerned with finite groups only. 
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(a) Time evolution of © with scalars in (b) Time evolution of © with scalars ex- 
Z2. tended to Z4. 



Figure 1: (a) is a projection of (b) induced by Z| -» l\. 

structure. We will then proceed, in section |3j to expose an algorithm 
taking as input the local transition rule and outputting a description of 
the spacetime diagram. This allows us to compute salient features of 
these fractals, such as their fractal dimension and their average color. 

1 DEFINITIONS 

1.1 Generalities on summable automata 

Monoids 

We want to discuss "summable automata", for which it makes sense to 
talk about the influence of a single cell on every other cell, and where 
the global transition function can be reconstructed by "summing" all 
these influences. So, if E denotes the alphabet, instead of the usual 
local transition function — » E, a summable automaton is naturally 
defined by a function E — >■ H 1 . What is then the minimal structure 
on E that would make such a definition work? These influences have 
to be "summed", so we need an operation on E. Since the strip is 
infinite, an infinitary operation would do, but that wouldn't give us 
much to work with. Instead, it seems reasonable to consider a binary 
operation +. In the same spirit, when we think of the superposition 
of influences coming from each cell, no notion of order between the 
cells is involved; even if in the one-dimensional case a natural order 
can be put on the cells, it would be less than clear what to do in higher 
dimensions. We require therefore that + be associative and commu- 
tative. The last requirement comes from the fact that, given only the 
global transition function, we want to be able to isolate the influence 
of one cell; that is why we demand that + have an identity element, 
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which makes now (E, +) an Abelian monoid. Of course, in order for 
all of this to be relevant, the transition function has to be a morphism. 

Let I be some finite subset of Z and / a morphism from E to E*. 
From / one can define the global transition function as an endomor- 
phism F of E z by 



/ E z -> E z 

r =( r »)nez ^ £/(>W)i 



(1) 



Let a be the right shift on Ez, i.e. c(r) n = r n _\. We have F o a = 
a o F, which means F is translation invariant. Also, F(r)„ depends 
only on the values r„_, for i E I; since I is finite, F is a one-dimensional 
cellular automaton on the alphabet E, with neighborhood included in 
— I. Conversely, if F is an endomorphism of E z defining a cellular 
automaton over the alphabet E, then one can choose a neighborhood 
I, and define, for i E I, 

f(s)i = Hs)u (2) 

where s is the word of E^ defined by s n = { f ^. , e denot- 

L e otherwise 

ing the neutral element of E. 

Groups 

We will now consider the case when E is a (finite abelian) group. For 
p prime, let E p be the subgroup of E of elements of order a power 

of p; then E is isomorphic to FJ Ep, and every endomorphism of E z 

V 

factorises into a product of endomorphisms of the E z 's. It is therefore 
enough to study the case of the (abelian) p-groups: let us assume E is 
a p-group. 

It is a well-known fact (see for instance section 1-8 of |7|) that E is 
isomorphic to x Z^ 2 x ■ ■ ■ x Z^ with > k^\ > . . . > k\ — 
k. Now let us consider an endomorphism a of E and let ej denote 
(0, . . . , 0, 1, 0, . . . , 0), where the 1 lies in position j. When i > j, there is 
a natural embedding s, y of Z^t ; into Z ^ , namely the multiplication by 

pkj-kj g g mce nas orc j er p k j / oc(ej)j E Z kj has to be in the image 
of su when i < j. We can therefore associate to a the endomorphism 
of Z^ given by the matrix A (a) E .^(Z^*) defined by A(a)j / j = 
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For instance, if G is Z32 x Z4 x Z2, and a is defined by a (1,0,0) = 
(3,2,1), a(0, 1,0) = (24,0,1) and a (0,0,1) = (16,2,0), then the corre- 
sponding matrix of ^#3 (Z32) would be 



A(«) = 



Let us give a summary of the construction we have just exposed. 

Proposition 1. For every finite abelian p-group G and endomorphism a. of 
G, there are positive integers k and d, an embedding s of G into Z^ t , and an 

endomorphism A(a) ofZ!^ k such that the following diagram commutes: 




(3) 



z* 4 J^U z* k 



This implies that to study the behaviour of CA on abelian groups, 
it is enough to study the case where these groups are of the form Z^. 

R-modules 

We will actually consider the more general case where R is a finite 
commutative ring, and E is a free R-module of dimension d, i.e. iso- 
morphic to R d . The first reason for doing so is that it does not compli- 
cate the mathematics. It will also appear more efficient to understand, 
for instance, F 2 4 as a 1-dimensional vector space over itself than as 
a 4-dimensional vector space over F 2 : the former simply bears more 
information, and therefore implies more restrictions on the form of a 
CA, so that more can be deduced. 

For any ring B, B [u, denotes the ring of Laurent polynomials 
over B; it is the ring of linear combinations of integer powers (neg- 
ative as well as nonnegative) of the unknown u. Applying this to 
B — Hom R (Z,), we can associate to the function / the Laurent poly- 
nomial r{f) E HomR (E) [u, u^ 1 ] defined by 

r(f) = E /(On""- (4) 

neZ 

t is an isomorphism of K-algebras between the linear cellular au- 
tomata on the alphabet E with internal composition rules (+,0) and 
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HomR (E) [u, , which can be identified with {K \u, ) be- 
cause E ~ R d ; we are going to think and work in this former algebra, 
so from now on a linear cellular automaton T = t(/) will be for us an 
element of (R [u, ) . 

1.2 Related work 

Many papers have been published about the fractal structure of cel- 
lular automata spacetime diagrams. We give here a short review and 
point out the differences to our approach. When we mention d and k 
we are referring to (Z k [u, w -1 ]). 

Il5l In this work, Willson considers the case d = 1, k = 2. In order 
to determine the fractal dimension of the spacetime diagram, he 
analyses how blocks of length n in the configuration of time step 
t are mapped to such blocks in steps 2f and 2t + 1, a technique 
we also use in section [3] 

1131 1141 Takahashi generalizes Willson's work to the case d = 1, with 
no restriction on the value of k. 

1 5 . 6 1 Haeseler, Peitgen and Skordev study the fractal time evolution 
of CAs with special scaling properties, the weakest of them be- 
ing "weakly p-Fermat", where p is some integer, which includes 
the case d = 1, k — p. Let us briefly introduce the p-Fermat 
property and show why the CAs that we study do not have to 
be p-Fermat. Let Tip be the scaling map 

A CA T is weakly p-Fermat if for all s 6 E, n 6 N and x 6 Z, 

T n P(s) x = e^TZ p T n (s) x = e. 

Let us now consider 

M - 1+ 1 1 + u )^ 2 (Z 2 [u,u^]). (6) 
We will use this example throughout the paper. It generates the 



time evolution depicted in Figure 2a A general nearest-neighbor 
p-Fermat CA produces a time evolution that reproduces itself af- 
ter p steps in at most three copies located at positions { — p; 0; p}. 
After 2p steps we have five copies at most. This creates areas 
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(a) Time evolution of © starting from £ = (b) Time evolution of a general nearest 
(J). neighbor p-Fermat CA. 



Figure 2: This figure shows that cannot be a p-Fermat CA. In a p- 
Fermat CA at least the white areas are filled by the neutral element e; 
has a different pattern. 



filled with the neutral element e shared by all p-Fermat CAs for 



a fixed p. In figures 2a and 2b we can easily see that does not 
exibit these areas; therefore it is not p-Fermat. Furthermore p- 
Fermat CAs that are not periodic are irreversible, while we also 
allow reversible CAs, being again one example. 

01 ED Allouche, Haeseler, Peitgen and Skordev study recurrences in 
the spacetime diagram of linear cellular automata, from the an- 
gle of A: -automatic sequences, which we will not define in this 
paper. However they require E to be an Abelian ring and the 
CA to be a ring homomorphism, which is again essentially the 
case d — 1. 



1101 ITU Moore studies CAs with an alphabet A on a staggered space- 
time, where every cell c is only influenced by two cells a and b of 
the last time step. The update rule is c = a • b. He requires (A, •) 
to be a quasigroup and studies different special cases. First let us 
note that these CAs are either irreversible or trivial, while ours 
don't have to be. Thus, although it is possible to bring our CAs in 
the form of staggered CAs, the results of Moore do not apply. In 
his setting, our CAs would be of the form c = a • b = f (a) + g(b) 
for some homomorphisms / and g. For (A, •) to be a quasigroup 
means 



Ma,b £ A 3\x,y G A a* x = b Ay • a — b. (7) 
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In our case, these equalities translate respectively as g(x) = b — 
f(a) and f(y) = b — g(a). The right-hand sides can be arbitrary 
elements of A, therefore / and g have to be isomorphisms, as 
indeed required in |10|. 

The angle of study of Moore is also different: he does not exactly 
study the fractal properties of the spacetime diagram, but rather 
the complexity of the prediction — "What will be the state of 
this cell after f steps?". Describing the spacetime diagram with a 
matrix substitution system is one way of proving that prediction 
is an easy task — for instance it makes it NC. 

1 8 1 Macfarlane uses Willson's approach and generalizes parts of it to 
some examples of matrix-valued CAs, including 0. However, 
the transition matrix is obtained heuristically — "by scrutiny of 
figure 9" — from the spacetime diagram, instead of being al- 
gorithmically derived from the transition rule (as in the present 
work). The conclusion (section 6) suggests that the analysis of 
is easily generalizable to matrices of various sizes over various 
rings, so in a sense the present article is but an elaboration of the 
concluding remark of QQ, although we have to say we do not 
find this generalization to be that obvious. 

The heart of our proof is in section|3] In a nutshell, whereas most of 
the techniques used in our article can be traced back to older articles, 
the new one that allows us to extend the analysis to a larger class of 



automata is the introduction of a in Equation 1 19 1. The idea in doing so 
is to get rid of the complicated noncommutative ring structure and go 
back to a simple linear recurrence, as stated in Proposition |4] Since a 
linear recurrence is precisely where the analysis started from, it could 
seem at first sight that nothing is gained in the process, but the new 
recurrence actually does not define a cellular automaton. Instead of 
defining line n + 1 from line n, it cuts right through to line mn, thus 
establishing a scaling property. 

1.3 Different CAs 

While we use 0, which has very special properties (being reversible, 
over a field of characteristic two, and described by a 2 x 2 matrix), as 
our example throughout the paper, the analysis applies of course to 
all other linear CAs. In this section we give a short overview over the 
variety of spacetime diagrams these CAs generate. Let us start with 
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small changes to 0. For our first example we keep m = k = 2, but 
change the determinant to u. We only change one entry of the matrix: 



©u 



+ 1 



(8) 



The spacetime diagram is displayed in figure 3a and shows how 
much difference a small change in the update rule can make for the 
spacetime diagram. 




(a) Time evolution of . (b) Time evolution of Tp 4 . 

Figure 3: Spacetime diagrams of non-clifford CAs. 



Let us now modify in a more subtle fashion: 



- 4 =(! «-i +!+«)• (9) 

The hidden difference with is the underlying ring, which has 
now been extended from Z2 to Z4. 0j.- = 4 contains in some sense more 
information than 0, since is induced from 0,t = 4 by the projection 
Z4 — » Z2. Consequently the spacetime diagram of is nothing but a 
projection of that of 0/ c= 4, as illustrated in figure [Tb| 

The last CA we want to present lies in ./#2 (^[u, w^ 1 ]), where F4 
is the finite field of order 4, here identified with F2 [a;] / {to 1 + co + 1 ) . 
The corresponding matrix is 

t F4 = ( A V do 



11 



(to + \)u 1 + CO + u 



If one wants to avoid calculations in F4, this CA can be translated to a 
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CA in (Z2 [u,u 1 ] ) , namely 



T, 
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u u" x + l 
1 1 + w / 



contains 



Its spacetime diagram, as can be seen shown in figure 3b 
patches of checkerboard pattern. Somehow, they trivialise most of the 
usual properties of the figure: for instance they make its fractal dimen- 
sion 2, even if the fractal structure can hardly be considered trivial. In 
order to access more interesting properties, it is possible to blank this 
pattern out, considering it as just "another shade of white". This can 
be trivially done on the matrix substitution system, by removing the 
states from which the blank state is inaccessible. 

1.4 Coloured spacetime diagrams 

The mainstream setting when studying the fractal structure of space- 
time diagrams is monochromatic; we introduce colors in the picture. 
Instead of considering simple compact subsets of the plane, we will 

have a finite set of colors c & and compact subsets of (IR ) . Let b £ c € 
be the additional "blank" color and c : E — > 'to U {b} a coloring of £ 
such that c(0) = b. To determine a colored spacetime diagram, we need 
furthermore to be given an automaton T £ (K [u, ), an initial 
state £ 6 R d , and an integer n. The corresponding colored spacetime 
diagram is then the rescaled diagram obtained by iteratively applying 
T n times on £. 

Formally, for n, i,j E N, let S n n be the full square centred in ~ (i,j) 
and whose edges, parallel to the axes, are of length ~. To each pos- 
itive integer n and color c £ ^ is associated a compact subset of the 
plane <^ n (c) which is the union of the S n n's such that < j < n and 
c (T1 — c. The colored spacetime diagram of order n is then the 
function &> n : c >—> £P n (c). A sequence of colored patterns (^ a n) nS ]N of 
spacetime diagrams is said to converge to some colored pattern if 
for every c E % ', (^ 2 n(c))„ e]N converges to ^ 2 o(c) for the Hausdorff 
distance. 

We can now state our main result. 

Theorem 1. Let G be a finite abelian p-group. For every cellular automaton 
over G that is also a group homomorphism, there exists a positive integer m 
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such that for every fixed initial state the coloured spacetime diagrams of order 
ptnn conver g e zvhen n goes to infinity. 

In general, to know about the fractal structure of a cellular automa- 
ton over some finite group G, write G as a product of p-groups and 
study each p-component of the spacetime diagram independently; ac- 
cording to Theorem [lj each component generates a fractal pattern. 
Then, since the logarithms of the prime numbers are rationally in- 
dependent, it is possible to find a sequence of resized spacetime dia- 
grams that converges towards a superposition of these different com- 
ponents with arbitrary independent rescaling coefficients, but there is 
no direct generalization of the theorem. For instance, even in the sim- 
ple case of Pascal's triangle modulo 6, there is no real number a. > 
such that the diagrams of order [ct n \ converge; however those of order 
t n will converge as soon as the fractional parts of log 3 (t n ) and log 2 (t n ) 
both converge, and then their limits determine the limit pattern. The 
situation is very briefly described in section 5 of [14J. 

Matrix substitution systems 

We will show how to find a suitable description of the limit pattern 
in the rest of this article. We now explain exactly what it means to 
generate a colored picture by rules of substitution, and how to take 
the limit of all these pictures. This is a generalization of the usual 
monochromatic description that can be found for instance in [9 15j|5l, 
and which corresponds in our setting to the case where all the colors 
are mapped to "black". 

Let V be a finite alphabet; because we want colors, compared with 
the usual definition of a matrix substitution system, we don't have to 
include a special "empty" letter. A matrix substitution system is then 
a function <3> : V — >■ y[ 1;r l ; for some integer r. Together with a set of 
colors c € and a coloring c : V — »■ c € , it defines colored patterns, much 
in the same way cellular automata do. With the previous notations, at 
each step n, the pattern is the union of squares S r n ( ; ; of different 
colors, for different i's and j's; each one of them is indexed by some 
letter in V. 

Then at step n + 1, each colored square of color c indexed by v £ V 
present in the nth step pattern is replaced by r 2 smaller squares that 
pave it; these smaller squares are given by 2#{v) and indexed accord- 
ingly. To such a matrix substitution system we can associate a multi- 
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graph T = ( V, E) where the set of vertices is V and we put as many 
edges from v to w as there are w's in Q){v). 

A plain matrix substitution system is one of the usual kind: no col- 
oring, and V contains a special letter e such that & (e) is a matrix full 
of e's and c(e) = b. In the multigraph associated to a plain matrix 
substitution system, e is excluded from the set of vertices. 

We want to generalize the usually property of convergence of the 
patterns defined by plain matrix substitution systems. This will be 
done by the conjunction of the two following propositions. Let us first 
remind some notions on graphs: the period of a graph is the greatest 
common divisor of the lengths of all the cycles in T; a graph is aperiodic 
if it has period 1. 

Proposition 2. If every strongly connected component of I is aperiodic, then 
(^")neN converges. 

Proof. To each color c e we associate the plain matrix substitution 
system obtained from 3) simply by turning some letters into £. For 
v G V, let X c (v) be the set of integers n such that there exists a path of 
length n in T connecting v to a letter of the color c. Since the strongly 
connected component containing v is aperiodic, X c (v) is either finite 
or cofinite. Those letters v 6 V such that X c (v) is finite are sent to e, 
and this defines @ c . If X c (v) is finite and v' can be reached from v, 
then X c (v') is also finite; therefore, 2l c is indeed a substitution system. 
Let M be such that for every v G V, either X c (v) or its complement is 
strictly bounded by M. 

Let us now compare two sequences of figures. The first one is 
(^ n (c)), the subpattern of color c defined by @. The second one is 
(^n)/ m e one obtained from £^ c ; we know that it converges to some 
compact By construction, ^ n+ u{ c ) is included in 3P c n , and for 
every black square of there is a black subsquare in £P„ + m{c). The 
Hausdorff distance between ^„ + m( c ) and therefore converges to 
0, so that (&>„ (c) ) converges to □ 

For a graph T, let T k = (V, E k ) where E k is the set of couples (v, w) 
such that there exists in T a path of length k from v to w. 

Proposition 3. For every (multi)graph T, there exists k such that every 
strongly connected component ofI k is aperiodic. 

Proof. Each strongly connected component A of T has a period p(A), 
so that AP( A ) is aperiodic. Let k be the least common divisor of the 
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p(A)'s; then each strongly connected component of T induces an ape- 
riodic graph in T °, but it is possible that, in the process, it broke down 
into several connected components, so that r*-° might not have the re- 
quired property. The procedure then has to be repeated from T k ° to 
obtain r^ 1 , and so on. Since the strongly connected components of 
P-o—fci+i are included in those of r* *' >, this process reaches a fixed 
point, which is a graph with the required property. □ 

Ergo, a colored matrix substitution system defines a convergent col- 
ored pattern when considering the steps that are a multiple of some 
well-chosen integer m. So, in order to prove Theorem [lj all we need 
to do is find such a substitution system. This will be done in a special 
case in the next section, and in the general case in section|3] 



2 A SPECIAL RECURSION SCHEME FOR 

The aim of this section is to give the most direct and natural expla- 
nation of the fractal structure generated by that we are aware of. 
Modulo some caveat, it applies effortlessly to all invertible elements T 
of Jii (R[u, m -1 ]), where R is a finite Abelian ring of characteristic 2. 
This section is not vital to the proof of the general case presented in|3J 
and can therefore be skipped by the impatient reader. 

We will deduce informally the basic structure of the spacetime di- 
agrams from a simple recursion relation for the 2"th powers of T. 
The characteristic polynomial of T, Pt(X), is equal to X 2 + (trT) + 
detT. According to Cayley-Hamilton theorem, Pt{T) = 0, so T 2 + 
(tr T) T + (det T) I = 0. Multiplying this equation by T , we get 
T = (det T) T -1 + (tr T) I. Let us denote T = (det T) T _1 , which we 
will name the dual of T; since we are in characteristic 2, by repeatedly 
taking the square of this equality, we obtain 

VneN T 2 " = f 2 " + (tr T) 2 " I. (11) 

Taking the trace of this equation, we get tr T 2 " = tr T 2 " ; in par- 
ticular, tr T = tr T so Equation jlTj is also valid when swapping T 
and T. Let It be a finite set, and the A,'s elements of R such that 
tr T = YLieh ^i 11 '- Then we have 

VneN (trT) 2 " = £ (A,) 2 V\ (12) 
iel T 
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We do not yet specify the initial state; as a matter of fact, it will 
prove to be largely irrelevant. The only thing we ask for now is that it 
is nontrivial (and finite). 

Consider for instance 0; we have det© = 1 and A, = X{-i ; 0;i}(0- 
We start with the spacetime diagram corresponding to 2" steps; it is 
rescaled to a triangle with vertex coordinates {(0,0), (—1,1), (1,1)}. 
Taking equations (TT) and (12) , we can see that the state at the 2"th 
time step can be decomposed into a sum of several copies of the initial 
state (the positions are governed by the coefficients of the trace) and a 
configuration that can be derived by applying T 2 " to the initial state. 
In the next 2" steps, this configuration will contract itself to the initial 
state, which is shifted according to log [( det T, as T is the inverse of 
T composed with the shift (det T)I. The copies of the original initial 
state evolve according to T. This is illustrated in Figure 13] The figure 




Figure 4: The whole figure is the sum of |Z| +2 parts. 



suggests to divide the spacetime diagram into four parts A, B, C, and 
D as shown in Figure 5a which overlap only on a single cell strip at 
the borders. 

Ai, A3, ad A4 are copies of A\, and V is marked by an upside down 
A because it is the reverse evolution of the initial state under the CA 
— so actually it should logically be named V or Y , but since it al- 
ways appears upside-down while A always appears straight on its 




(a) The first substitution rule. (b) The second substitution rule. 

Figure 5: The first and the second substitution rules. 
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feet, there is no risk of confusion. 

Let us assume that the sequence of rescaled spacetime diagrams 
up to step 2" actually converges. Then that means A\ should be, in 
the limit, a copy of the whole picture A, downsized by a factor 2, so 
we rename it A. This gives us the first substitution rule, represented 



in Figure 5a The other three parts are still unknown, and we will 



name these patterns B, C and D. Equation 11 tells us what the other 
substitution rules are. 

Since Equation [ll] remains true after swapping T and T, V admits 
likewise a partition into V, Q, D and Q. Summing all the parts shown 
in Figure |4j we get the top rightmost pattern of Figure [9] Superimpo- 




Figure 6: Second step of the decomposition. 



sing our first substitution rule (Figure 5a I with our second step of the 
decomposition (Figure pi, we get the new substitution rules repre- 



sented by Figures 5b 7ajand 7b 



All the other substitution rules are deducible from these ones. First 
they are linear: for instance the substitutions for C + D is the sum of 



the substitutions for C and D, as shown in Figure 8a We don't know 
a priori what the sum of two patterns is, but we know that summing a 
pattern with itself should give 0, for we are working in characteristic 




(a) The third substitution rule. (b) The fourth substitution rule. 

Figure 7: The third and the fourth substitution rules. 
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2. In this case A + A and B + B cancel out. In figure [9] one can see 
how the characteristic white spaces emerge from the above substitu- 
tion rule. 




(a) The fith substitution rule. (b) The sixth substitution rule. 



Figure 8: The fifth and sixth substitution rules are the superpositions 
of the third and fourth and the first and second rules respectively. One 
can see, how white spaces emerge because even superpostitions of the 
same pattern cancel out. 

The problem with this scheme is that what is happening goes be- 
yond simple juxtaposition of patterns. There can be cancellation at the 
border between patterns. And, sure, we know C + C is blank, but how 
do we know C + D, for instance, is not? Well, generally we don't. If 
the initial state is itself blank, then the whole figure would be, and all 
tiles being blank is certainly a fixed point for all the substitution rules. 

In this case, however, everything turns out well. It should first be 
noticed that in every part of the picture not tagged as "blank", an 
A pattern can be found by refining a few more steps. Formally, let 
G = (V, E) be the graph whose vertices are the different possible tiles 
(i.e., in this case, A, B, C, D, V, Q, D, Q, and the sums modulo 2 of 
tiles having compatible shapes, including the blank tile 0), and edges 
represent the transition rule in the following way: each vertex has four 
edges coming out of it, each one pointing to one of its subfiles. In our 
case, the graph has the property that the set of vertices accessible from 
A, minus 0, form a strongly connected component. 

We may then distinguish two cases: either A has a point in its in- 
terior, or it has points only on its border triangle. In the first case, the 
unique non-empty compact defined by the substitution rule is actu- 
ally the figure we're looking for. Indeed, it follows from the property 
of connexity — cf. Proposition [2] — that every non-zero tile actually 
appearing in the decomposition of the figure has a non-empty inte- 
rior. Thus, no matter what happens at the boundary between tiles, the 
figure constructed this way will always converge to the same compact. 
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Figure 9: Three decomposition steps for the spacetime diagram of 0. 
One can clearly see the characteristic white spaces emerging. 



3 RECURSION AND MATRIX SUBSTITUTION SYSTEM 

We will now present a general method to calculate the fractal dimen- 
sion and average colour of the spacetime diagrams of linear CAs in 
[u, m -1 ]). We will again illustrate the method with our exam- 
ple 0, whereas the derivation is carried out for the general case. Thus 
the algorithm works as well on all CAs obeying our definition, e.g. the 



CAs presented in section 1.3 as it works on 0. Of course with larger 
neighborhoods and groups of higher order the substitution system be- 
comes larger and larger, so that one might want to use a computer to 
derive the substitution system. 

Our approach is the following: from the minimal polynomial TI of 
the CA T (or any other polynomial fulfilling n(T) = 0) we derive a 
recursion relation for the Tf 's, the coefficients in u x of T y . We then 
forget about every other piece of information we might have on T, 
to concentrate only on this recursion: this shows that the fractal struc- 
ture, except for contingent blank spaces, can be essentially derived just 
from the minimal polynomial of T. We further develop our recursion 
scheme for T until we can express every T\ in terms of the r, of the 
first m time steps with coefficients a.; (x — i,y). With a simple grouping 
of cells we deduce a matrix substitution system that enables us to gen- 
erate the spacetime diagram of step t = k" +1 directly from step t = k n . 
Using this substitution system we can calculate the fractal dimension, 
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the average colouring, and given an initial state also the whole space 
time diagram. 

Now let n(X) £ R[u, [X] be a monic polynomial such that 
n(T)=0: 

m— 1 

n(X) = X m - £ An,/* 7 '- (13) 

/=o 

According to the Cayley-Hamilton theorem, which we can apply in 
our case because T is an endomorphism of a finite-dimensional free 
module over an abelian ring (see Theorem 3.1. of [7]), the characteris- 
tic polynomial of T fulfills this condition, therefore we can always find 
such a polynomial. Let I be the finite set of exponents z's such that the 
coefficient in u 1 of n, seen as an element of R[X] [u, m -1 ], is nonzero, so 

that we can write An , = ; i u '- ^ is n °t to be confused with the 

ieX 

neighbourhood of the CA, I, which won't play any role from now on. 

For any x,y,n £ N, (x + y)P = x p + y p [p], and if x = y[p"] then 
x? = yf[p n+1 ]. Therefore, for any x,y £ R and n £ N, (x + y) p + = 
[xP n + yP n y l 1 . Since the powers of T commute pairwise, we get 

„n+2(/-l) m ,p"+'- 1 W.+H,y 



For each z,/, the sequence (A^ . .) is ultimately periodic. There exist 

„M+N „M 

therefore integers N and M such that for all Aj-^ ■ = Aj-j- •. Let 
k = p N ; substituting n by M + Nn in ( 14 1, we get 



Hence, if we note m' = p M+2 (' _1 )m and expand this equation, we 
find that there is some finite subset I' of Z and some elements ^ of 
R, for z £ X 1 and £ [[0; m' - 1], such that for all n £ N, 



m'-l 



7* ,m '= £ 1] m,ju m T kn i. (16) 



18 



We have now used everything we needed to know from the multi- 
plicative structure on the ring of matrices. As announced at the end 
of section [L2| we will now get rid of it and concentrate only on the 
linear recurrence relation that we have just derived. Remember that 
V £ JK d (R [u,u- 1 ]), and we are interested in the coefficient of Tl in 

iij, denoted so that T> = £ V-u 1 . We thus get the following rela- 

ieZ 

tion: t\ m +J/ = YjieX' UlLo 1 ^hj^x-k"!' wn i cn we rewrite in this form: 

Tt= E HjTf^ (17) 

where fu(y) = —k"i and gu(y) = y — k n (m' — j), which of course 
works with any n, but we will choose n = [logj. j^J . In order to em- 
phasise that the rest of the proof will use only a minimal structure, we 
state in the next proposition what will actually be proven, and change 
the notation from T, which was an element of (R[u, m^ 1 ]), to E, an 
element of a more arbitrary R-module. It is straightforward to check 
that T fulfills the hypotheses of the proposition. 

Proposition 4. Let M be a finite R-module, k a positive integer, A a finite set 
of indices, and for i £ A, }ij 6 _R, ft : fin; +oo[[— ^ Z and gj : fin; +oo[[— » N 
such that for all y G \m; +°o[ and t £ JO; k — 1], 

• gdy) < y; 

• fi(ky + t)= kfi(y) and gi (ky + t)= kgi{y) + t. 
For x £ Z x N, let 6 M be such that when y >m, 

* = Li<L- < 18 > 

Then there exists a finite set E and a function e : Z x N — > E such that 

• is a function of e(x,y); 

• for s,t £ JO; k — 1], e (kx + s, ky + t ) is a function ofs, t, and e(x, y). 

The introduction of a new function e in this proposition comes from 
the need of a scaling property, expressing that the state at point (kx + 
s,ky + t) can be deduced from the state at point (x,y). Such a prop- 
erty does not follow immediately from Equation jTH) , but it is possible 
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to expand the state space from M to E, and to put more information 
into e than into S, so as to fulfill the scaling property. An immediate 
consequence of this proposition is that the spacetime diagrams of 
of order k n can be described by colored matrix substitution systems, 
so that Theorem [T] will follow from Proposition [3] Let us now prove 
Proposition |4] 



(1,3) 




(0,2) (1,2) (2,2) (1,1) 

(-1,1) (0,1) (1,1) (0,0) (0,1) (1,1) (2,1) (1,0) (1,1) (2,1) (3,1) (2,0) 



Figure 10: Recursive calls to ( 18 for T = 0. 



If y > m, we can recursively apply Equation < ji~8) to give an expres- 
sion of S^- in terms of a linear combination of E^,'s with y' < y. 

Let us consider the case T = as an example. In this case, R = F2, 
k = m = 2, and FI (X) = X 2 + (w -1 + 1 + u) X + 1. Equation (181 



gives E^ = Sj + Eg + E 2 + S^. Whereas the decomposition halts here 
for Si, it goes on for the other three terms of the sum, again with n = 0. 
For instance, Eq is transformed into Sjj + E^ + Eg + S]. This gives 



the tree of recursive calls in Figure 10 where each node is labelled 
(x,y) for E^; in that particular case, because R = F2, the only possible 
coefficient are and 1, so they are simply represented by the absence 
or the presence of the corresponding term. 

Starting from any point (x, y) eZxN and using recursively equa- 
tion (18), we get to the expression 



m— 1 



E E a n,;,/(*,y)S-, 

i ;=0 



(19) 



10 



this corre- 



which we take as a definition of arj,j ;{x,y). On Figure 
sponds to (parity) counting the leaves tagged (/,/); for instance, there 
are 4 leaves tagged (1,1), so an ,i,i(l,3) = and (1,1) does not ap- 



pear as a node in Figure 10 



Since the relation | [18) is invariant under translations of the param- 
eter x, we have otu(x, y) = ccq^x — i,y). Noting otj := ocq^, we get the 
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(1,3) 



(0,0) (1,0) (2,0) (-1,1) (3,1) 



Figure 11: E^ as a sum of the simplest terms, when T = &. 



following equation: 

m— 1 

s? = EE */(*-*/y)sj. (20) 

iez j=o 

Let us now show that a. (kx + s,fcy + £), for s,f G [0;fc — 1], is a 
function of s, t, and the a. (x',i/)'s, where x' ranges over some neigh- 
bourhood of x. By substituting x with kx + s and y with fcy + t in 
Equation (18) , we get 

^kx+s - }*i'-'1cl x+ f i (y)) + s- 
i'6A 

Because the indices and exponents of E on the left and right side 
of this equation have undergone the same transformation (x,y) i— >■ 
(kx + s,ky + t), we arrive recursively at this point 

m— l 

EE «/(*- ^s:; ™ 

z€Z ;=0 

which we want to compare to the following equation, directly de- 
duced from ((20}: 

m-l 

= E E ^ + s-i,ky + t)E\ (23) 

iez j=o 



Of course, there can be terms in p2| with kj + t > m, so that the 
decomposition is not over: it then needs to be performed to its end. 
For instance, in our example T = 0, if instead of the tree of recursive 
calls for E^ we want that of Ej = Sjjxl+l' we ' <now ma * me upper part 



of the tree will be Figure 12 



As for the lower part, the decomposition is not over, so it needs to be 
performed to its end. Thanks to linearity, we don't need to remember 
the whole tree, just its leaves. As is well known, a tree grows from its 
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(3,6) 




(1,4) (3,4) (5,4) (3,2) 

(-1,2) (1,2) (3,2) (1,0) (1,2) (3,2) (5,2) (3,0) (3,2) (5,2) (7,2) (5,0) 



Figure 12: Recursive calls to (18 for T = &, minus some leaves. 



leaves, and in the same way we need only apply the transformation 
(i,f) i — y (2i + 1, 2j) to the tree in Figure 11 in order to find the correct 
decomposition for S^, which gives Figure 
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(-1,2) (7,2) 

(-2,1) (-u) (oa) (-y)) (64T^PSiT~ro)) 



Figure 13: Trees grow from their leaves. 

However, even keeping track only of the leaves would require an 
unbounded alphabet, that is why we defined a.n j(x, y) = ocjijjix + 
i,y), i.e. from the tree for we want to keep track only of the leaves 



labelled by (0,/). Alas, that will not work so easily: in Figure 13 (0, 1) 
appears once as a son of (—1,2), which comes from node ( — 1,1) in 
So it appears we cannot just forget about everything else. 
As we just saw on one example, an,- (kx + s,ky + f) typically does 



Figure 



not depend on an,- (*/]/)• Instead the final decomposition of 1 22 will 
relate the coefficients &nj(kx + s — i,ky + t) of the E^ in 1 23 1 to sums 
of the oc^jj(x — i,y)'s. Therefore au, (kx + s,ky + f) depends on the 
an,- (x + i, y)'s for i ranging over some finite set 5? depending only on 
n. However, this is not much of a problem, as a simple grouping will 
take care of it — a technique commonly attributed to ||T5] . Let us show 
that 5? is finite. Since / £ {0; . . . ; m - 1} and t G {0;...;k-l}, the 
kj + t appearing as an exponent of E in |22| is in {0; . . . ; km — 1 } . We 
therefore have to use at most (k — 1 ) m recursive calls to (18) in order 



to get down to coefficients E^ with y < m, each one of them decreasing 
the exponent by at least 1 . Each one of them also increases the index by 
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fi(y); since both A and 10; km — 1] are finite, the set of possible //(j/)'s 
is also bounded by some M, and the total variation in the index, i.e. 
y, is then bounded by (k — 1) mM; let US Say y C | timing rfmax] • 
Let us now introduce ^n r (x,y) = (a.(x — z,y)) !(Ei y/, where 5?' = 



Pmin; ^max]/ Vax being such that <5 n 



and S m ;„ such 



that Vin < 



This time, for s, t 6 {0; . . .;k — 1}, 



j6n,. (kx + s,ky + t) does really depend only on /5n,. (*/ J/)- Indeed, 



j6 n ,- (be + s,fcy + f) = (a. (be + s — i,ky + 0)igj" / 



and each a. (kx + s — i,ky + t) depends on (a. I v 
only; the choice of has been made so that j 



S — 1 

k 



s—i 
k 



eludes the proof of Proposition |i] as we can choose E = Ml° ; 



G S". Thiscon- 



with e(x,y )(;',/) 



i,y). 



3.1 Example: O 

In the case of 0, Equation p2} becomes 



(24) 



The first term is now elementary, but the second one has to be de- 
composed once more, i.e. 



(25) 



which is the end of it if t = 0, but not if t — 1, where we get 



7?3 _ 77O . 7^0 ,770 ,771 , 77I 

"2i'+s — "2i+s-l 1" ^2i+s 1" ^2i+s+l "r "2/+S-2 "r l -'2i+s+2- 



(26) 



Comparing this with 1 23 1 we can deduce the substitution rule of a@ . . 
It can then be written in the following way, where for convenience a© 



is represented as 



«©,1 
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a ,.(2x,2y + l) 


a 0/ .(2x + l,2y + l) 


a .(2x,2y) 


«e,.(2x + l,2y) 



II 



«©,o (*,y) + *e,i - 1/3/) + *e,i + 1/3/) 
a©,i(x,y) 




a©,i(^y) + + !'!/) 


a©,i(x,y) 
«e,o (x,y) + « eA (x,y) 


«©,i( x '3/) + a e,i(* + 1 >y) 




If we follow exactly what has been said in the general case, we 
ought to consider the grouping [—2; 3]. However, this general bound 
is obviously too rough in the case of 0, where we will just have to take 
[— 1;2]. We will represent the grouping in the form 



a 0/1 (x-l,y) oi @iX (x,y) a 01 (x + l,y) a 01 (x + 2,y) 
a 0/O (x-l,y) fl 0/O (x,y) g 0O (x + l,y) a: 0O (x + 2,y) 

The alphabet has thus size 256, and the substitution system is de- 
scribed by 



abed 
e f g h 






a+c+f b+d+g 
a + b b b + c c 


a+c+f b+d+g 
b b+c c c+d 


a+b b b+c c 
b+f c+g 


b b+c c c+d 
b+f c+g 



Let us denote by A the matrix having a 1 in position a and else- 
where, B the matrix having a 1 only in position b, and so on. For these 
matrices, we will denote the sum of matrices by a simple juxtaposition: 
AB will mean A + B, as the matrix multiplication has no meaning in 
this context. 



Since = 8 Xi qTq, the starting position, with which we describe the 



whole line number 0, is 





H 


G 


F 


E 





• • • . Since we 


have, for instance, the rule F — > 


B 


A 


-, the graph derived from 


F 


E 



this substitution system is aperiodic; that means that, in whatever way 
A, B, C, . . . are represented, either as colored dots or as white dots, the 
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pattern converges, and the fractal structure is described by this matrix 
substitution system (see Section |i~4) . 

To calculate the fractal dimension of our spacetime diagram we use 
the transition matrix of the matrix substitution system, which con- 
tains the information about the images of all states. The line corre- 
sponding to F would contain a 1 in the rows of A, B, E, and F and 
zeros elsewhere. As every cell gives rise to four new cells the sum 
of all entries in each column of the matrix is 4. We thus deal with a 
sparse 256 x 256 matrix. The base 2 logarithm of the second largest 
eigenvalue of this matrix is the fractal dimension of the spacetime di- 
agram (cf. for instance l!T5l ). Here this gives a fractal dimension of 
log 2 - 1-8325, as also found in 0. 

Let us note that up to this point our analysis for is word for 
word valid for all CA in ^#2(^2[ M / U of determinant 1 and trace 
w -1 + 1 + u. The additional information is only used for the actual 
colouring of the picture. In general all CA with the same minimal 
polynomial have the same substitution system, and in dimension 2 the 
minimal polynomial is entirely determined by the trace and the deter- 
minant. The fractal we get if we use the substitution system starting 



from ■ ■ ■ 





H 


G 


F 


E 





• • • is shown in Figure 


14 



In the case of the connection between the substitution system and 
the colored picture is very simple; let us take £ = L) as the initial state. 
Then the state of cell x after y iterations is 

®li = E(«0,o(^-^y)©? + «0,i(^-iy)©9^ 

' - (27) 

= a 0/O (x,y)£ + y> 0/1 (x-z,y)0i£. 

i 

Since ©I = ( J ! )' 1 = -i = ( o J ) and ©/ = when 

i £ { — 1;0; 1}, we have if £ = ( ^j*'^ ). This gives us a color 

V a 0,i (*>!/) / 

assignment for each state of the matrix substitution system, which cor- 
responds to simply dropping all states that include neither B nor F. 

We can now determine the average hue of the spacetime diagram 
making use of the eigenvector corresponding to the second largest 
eigenvalue of the transition matrix |l5l. Let us say (J), (J) and (\) 
are respectively coded by the colors c\q, cq\ and cw, let Cq be the white 
color. We determine which symbols of the alphabet belong to each 
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Figure 14: The general fractal time evolution of a linear CA with 
k = m = 2, determinant 1 and trace u^ 1 + 1 + u. Only the areas where 
the whole group of as is is marked white. Thus the image appears to 
have less white than the colored picture. In the limit of infinite recur- 
sion this effect vanishes, thus the fractal that is generated is actually 
the same. 



of the colors by looking only at the part (* ' o ^^|)- Then we just add 
up all the weights of symbols with the same color in the eigenvector. 
We get the following unnormalised coefficients: c 10 : 2(4 + a/17), c m : 
2(4 + y/Y7), and c n : 5 + y/Y7. 



In figure la this color code was used: Cio = , Cqi = H and c\\ 



I. We must therefore have the following average hue: 



CONCLUSION 

We have shown that every cellular automaton inducing a morphism 
of abelian groups produces a spacetime diagrams that converge to 
selfsimilar patterns that are superpositions of elementary patterns ob- 
tained for groups of prime power order. We exhibited an algorithm 
taking as input the local transition rule and outputting a description 
of these elementary patterns. We only studied the one-dimensional 
case in this article, but the analysis can be carried over to higher di- 
mensions with not much more ado. Instead of (R[u, w -1 ]), a n- 
dimensional linear cellular automaton would then be an element of 
^^(R[u-\, u^ 1 , «2, u^ 1 , • • ■ u n , m" 1 ]), matrix substitution systems would 
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become (n + 1) -dimensional array substitution systems, the system of 
indices in section 3 would be further complicated, and the spacetime 
diagrams would be harder to display. However, the generalization 
does not present any theoretical difficulty and makes nice pictures. 



For instance, Figure 15 is a view of the the fractal structure pro- 
duced by Q = ( ? ,„ , i . 

1 1 (1 + « + «- 1 )(l + i? + z;- 1 ) 



a variation on 



belonging to .Mi (Z2[w, u 



,v,v 



For obscure historical reasons, 



pyramids used to be constructed with their points facing up, so to con- 
tinue the tradition, this fractal is represented upside-down compared 
to the two-dimensional figures. Its neighborhood being a square cen- 
tered on (0,0), the pyramid has unsurprisingly a square base. We cut 
it open so as to reveal its insides. The sector that is cut out is delimited 
by two half -lines, the one that is visible on the right along (1,0), the 
sector having an arbitrary angle of approximately | n. 




Figure 15: A modern funeral labyrinth. 

For any positive integer n, the coefficient in v n of O", seen as an 





element of (Z? W, u x l) \v,v 1 1, is ( „ " 1x „ I.There- 

V L IJ 1 ' \ (1 + U + M" 1 )" J 

fore, the pattern on the face of the pyramids is given by (1 + u + 
« _1 ) 6 Jt\ (Z2 [«, u^ 1 ] ) , as can be indeed noticed on the figure, where 
half of this symmetric pattern is clearly visible on the external face. By 
contrast, the coefficient of Q" in v° is exactly 0", which explains why 



the cut along (1, 0) in Figure 15 looks like the half of Figure 14 (turned 
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upside-down). 

Here are some open questions and possible future developments: 

• Presumably, the m in Theorem [T] can always be taken to be 1. 
This is known to be true in the cyclic case, i.e. when d = 1, cf 

• The algorithm presented in this article, producing a description 
of the spacetime diagram in the form of a matrix substition sys- 
tem, has a high complexity, due to the large size of its output. Is 
this a necessary evil, or can more efficient descriptions be found? 
Could for instance the more elegant triangle-based substitution 
scheme presented in section |2]be naturally generalized? 

• To what extent can the algebraic structure be weakened? In- 
stead of the alphabet being an Abelian group, could we consider 
an Abelian monoid? Is it possible to get rid of commutativity 
and / or associativity? 
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